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Abstract. Disk galaxies with a spheroidal component are known to host Supermassive 
Black Holes (SMBHs) in their center. Unequal-mass galaxy mergers have been rarely stud- 
ied despite the fact that they are the large majority of merging events by number and they 
are associated with the typical targets of gravitational wave experiments such as LISA. We 
perform N-body/SPH simulations of disk galaxy mergers with mass ratios 1:4 and 1:10 
at redshifts z=0 and z=3. They have the highest resolution achieved so far for merging 
galaxies, and include star formation and supernova feedback. Gas dissipation is found to be 
necessary for the pairing of SMBHs in these minor mergers. Still, 1:10 mergers with gas 
allow an efficient pairing only at high z when orbital times are short enough compared to 
the Hubble time. 
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1. Introduction 

Observations of nearby galaxies show that 
Supermassive Black Holes (SMBHs) rang- 
ing between ~ 10* and ~ 10'' inhabit 
the centers of virtually all massive galactic 
spheroids, from massive ellipticals to pseudo- 
bulges of late-type galaxies. Their masses, 
as infeiTed from dynamical measurements, 
appear to be correlated with various properties 
of their hosts, e.g. bulge lu mino si ty an d 
mass (e.g. [K ormendv & Richs tond (Il995 )' 
iMagorri an et all (il998i) : iMarconi & Hunt 
(^2003)), velocity dispersion (iTrem aine et aL 
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I2OO2L iFerrarese & Merritd l2000h. c oncentra- 
tion of the light profile dGraharn et al. 2001). 
In the ACDM model, galaxies build up their 
masses hierarchically starting from initial, 
gravitationally amplified fluctuations (e.g. 
IWhite & ReesI (Il978h '): therefore, every time 
two galaxies merge, the remnant is expected 
to host two (or more) SMBHs. The formation 
of a binary SMBH sistem has been shown to 
proceed quickly once both the compact objects 
are embed ded in a circumnuclear gaseous 
disk (see M aver et alj (2007); see also Mayer, 
Kazantzidis & Escala in this volume), but 
whether the large scale merger can lead them 
to such a favorable configuration is still a mat- 
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ter of debate: while mergers between galaxies 
of equal mass have been quite extensively 
explored in literature and se em to lead to the 
forma t ion of a SMBH pair dKazantzidis et all 
l2005t ISpringel et all l2005h . little attention 
has been paid to the fate of SMBHs in events 
which are much more frequent in the typical 
history of a ACDM galaxy, i.e. minor mergers 
with mass ratios fr om 1:4 to 1:10 and less 
dStewartet al.ll2008l) . 

2. Simulations and results 

2.1. Initial Conditions 

We discuss here results from N-body/SPH 
simulations of mergers between disk galax- 
ies with mass ratios q - 0.25 and 0.1 at an 
unprecedented resolution. The galaxy models 
were initialized as three-component systems 
comprising a Hernquist bulge, an exponential 
disk with a gaseous mass fraction fg, and a 
dark matter halo with an adiabatically con- 
tracted NFW profile, plus a central collionless 
particle representing the SMBH, whose mass 
was chosen following the M rh - Mbuige rela- 
tion dMagorrian etalJll998h . For more details 
on the setup of the initial conditions for th e 
reference model see Kazantzidis et al ] d2005h . 
We also ran mergers with initial c onditions 
rescaled to z = 3 according to the 'Mo et_alj 
dl998h (MMW) and Bullock et al. (2001) mod- 
els, and keeping Vvii fixed, as expected for 
the high-z progenitors of o u r present-day m od- 
els dWechsler et al.l 120021: iLi et al.1 l2007h . A 
large fraction of the gravitational wave sig- 
nal from coalescences of SMBH binaries is 
expected to come from this cosmic epoch 
dSesana et al.ll2005l) . if SMBH pairing is effi- 
cient. The satellite galaxies are initialized self- 
consistently with the same three-component 
structure of the primary, and a mass in each 
component scaled down by q. The baricen- 
ters of the two galaxies in each run were 
initally placed at a separation equal to the 
sum of their virial radii, on parabolic orbits 
with pericentric distances equal to 0.2 times 
the virial radius of the more massive halo 
dKhochfar & Burkertll2006l: Ib enso nl l2005l) . We 
ran collisionless ("dry", with fg = 0) and gas- 



dynamical ("wef", fg = 0.1 and 0.3) simula - 
tions using GASOLINE dWadsley et al.ll2004 . 
We will define two SMBHs a "pair" if their rel- 
ative orbit shrinks to a separation equal to our 
force soften i ng (~ 20 pc); at these distances, 
iMaver et al.l (|2007|) have shown that sinking 
proceeds very quickly and a SMBH binary can 
be formed in ~ 1 Myr Since our main purpose 
is to check which processes favour or inhibit 
this pairing, we do not employ any additional 
prescription in order to to keep them at the cen- 
ter of their galaxies' potential wells or to facil- 
itate their orbital sinking. 

2.2. Collisionless Mergers 

The galaxy merger and the sinking of the 
lighter SMBH towards the more massive one 
can be roughly divided in three stages: dur- 
ing the first, the orbit of the satellite decays 
because of dynamical friction on the halo of 
the primary; the second encompasses the most 
crucial phase of mass stripping in the denstest, 
baryon-dominated region where the fate of the 
SMBH pair is decided; in the final stage the 
remnant settles to its final dynamical state. 

For q - 0.25, dynamical friction on the 
dark matter halo of the primary is efficient, 
and the satellite galaxy sinks down to a few 
~ 10 kpc from the center. However, a colli- 
sionless system is not able to dissipate the en- 
ergy gained through tidal shocks at pericen- 
tric p assages (iGnedin et al.lll999l:lTaffoni et al.1 
I2003h : thus, the satellite is disrupted before its 
SMBH can reach the center of the remnant. 
No pair is formed, and the smaller SMBH is 
left several kiloparsec away from the other; at 
these distances its dyn amical friction timescale 
dChandrasekharl[T943h is longer than a Hubble 
time, but this is also where gas dynamics can 
affect the sinking (see l2.3l l. 

On the other hand, even after 10 Gyr the 
1:10, z = dry merger is not able to enter the 
regime where the inclusion of gas could speed 
up the sinking, because dynamical friction in 
the halo is too slow. Mergers at z = 3, instead, 
allow for an efficient sinking on much shorter 
timescales. The MMW scalings predict, for a 
given Vvii, masses and radii that are a factor 
of Hiz - 3)/Ho ~ 1/5 smaller; as a conse- 
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Fig. 1. Color-coded density maps for stars (upper row) and gas (lower row) at the second apocenter for 
three different ^ = 0. 1 mergers at z = 3: from left to right, the case with gas gooling (/j = 10%), star 
formation (J^ = 10%) and star formation (/j = 30%) are shown. Each frame is 30 kpc across. In the runs 
with star formation, supernova feedback creates a clumpy and diffuse ISM, resulting in strong ram pressure 
stripping when the orbit of the satellite passes through the gaseous disk of the primary galaxy. 



quence, orbital periods associated with orbits 
of same energy and pericenter (in units of 7?vir) 
are reduced by the same factor. Therefore, the 
1:10 collisionless merger at z = 3 is completed 
in ~ 2.5 Gyr Like for q - 0.25, a wandering 
SMBH is left at ~ lOkpc from the center of the 
remnant. 



2.3. Gas Dynamics 

Gas dynamics can change the orbital evolution 
only during the second stage of the merger, 
when the satellite is moving through the dens- 
est, baryon-dominated region of the parent 
galaxy; this is when tidal shocks and ram pres- 
sure become so strong that gaseous dissipation 
can significantly affect mass stripping around 
the smaller SMBH. 

For this reason, the dry and wet, q = 0.25 
mergers differ only after the first couple of or- 



bits (t ~ 6 Gyr) . As a lready pointed out in 



iKazantzidis et al. (l2005h . gaseous dissipation 
is a key element in the second stage of SMBH 
pairing in unequal-mass mergers: for 1:4 merg- 
ers, the presence of gas is a necessary and suf- 
ficient condition for the formation of a pair at 
the force resolution. The pair is embedded in a 
nuclear gaseous disk of radius ~ 500 pc, sup- 
ported mainly by rotation (~ 200 km s"'), and 
comprising a mass of ~ 3 • lO** Mq. 

In the 1:10, z - 3 merger with cooling, the 
gas in the outer disk of the satellite gets tidally 
stripped; however, most of it, after the first cou- 
ple of pericentric passages, is funnelled by tidal 
torques to the center of the satellite, ending up 
in a ~ 7 ■ 10^ M© circumnuclear structure. This 
central overdensity steepens the mass profile 
of the satellite, allowing it to survive subse- 
quent tidal shocks until it gets dragged down 
to the nucleus of the remnant, where a SMBH 
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pair is formed. The gaseous nuclear structure 
in the remnant is disk-Uke, of thickness equal 
to our softening length and supported by rota- 
tion in the disk plane (with rotational velocity 
~ 180kms-i). 

2.4. Star Formation 

We discuss now preliminary results from the 
last set of simulations, which include star for- 
mation and feedback from supernovae accord- 
ing to the prescriptions and fiducial parameters 
detailed in lStinson et al.i ( .2006.) . A comparative 
view of three q = Q.l, z = 3 mergers is shown 
in figure [T] 

When star formation is included, the 
Interstellar Medium (ISM) in the disks shows a 
multiphase and irregular structure, with larger 
scalelengths both in the radial and vertical di- 
rection. For this reason, the first pericentric 
passage of the satellite is not able to excite 
strong, coherent inflows which can steepen the 
density profile. Thus, during the second orbit, 
ram pressure exerted by the ISM of the primary 
galaxy strips most of the gas away from the 
satellite. The fraction of gas remaining in its 
center can then either be consumed by star for- 
mation, or stripped when the orbit reaches the 
densest regions of the primary: in this case, the 
small SMBH might not be able to pair with the 
more massive one before tidal shocks disrupt 
the core of the satellite. 

3. Conclusions 

Understanding the formation of SMBH bi- 
naries is of fundamental importance for the 
search of gravitational waves as well as for 
all studies of black hole demography and host 
galaxy coevolution. SMBH pairing in unequal- 
mass mergers depends crucially on gasdynam- 
ical eff'ect: satellites are disrupted too quickly, 
leaving wandering SMBHs in all the colli- 
sionless cases we studied. The presence of 
gas seems to be necessary for the pairing of 
SMBHs. While it also appears sufficient for 
q - 0.25, lower mass objects are more heav- 
ily affected by feedback from star formation, 
and their outcome could be very sensitive to 
the details of the physical processes involved. 
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